Electromagnetic radiation by gravitating bodies 

Iwo Bialynicki-Birulqj 

Center for Theoretical Physics, Polish Academy of Sciences 

Al. Lotnikow 32/46, 02-668 Warsaw, Poland and 

Institute of Theoretical Physics, University of Warsaw, Warsaw, Poland 



oo 

o 

O 

>>: 



Zofia Bialynicka-Birula 

Institute of Physics, Polish Academy of Sciences 

Al. Lotnikow 32/46, 02-668 Warsaw, Poland 

Gravitating bodies in motion, regardless of their constitution, always produce electromagnetic 
radiation in the form of photon pairs. This phenomenon is an analog of the radiation caused by 
the motion of dielectric (or magnetic) bodies. It is a member of a wide class of phenomena named 
dynamical Casimir effects, and it may be viewed as the squeezing of the electromagnetic vacuum. 
Production of photon pairs is a purely quantum-mechanical effect. Unfortunately, as we show, the 
emitted radiation is extremely weak as compared to radiation produced by other mechanisms. 
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I. INTRODUCTION 

The aim of this paper is to calculate the probability 
of photon emission by gravitating orbiting bodies. There 
are two separate mechanisms responsible for this process: 
electromagnetic and gravitational. The electromagnetic 
mechanism has its source in the variation of the material 
coefficients with time. Such changes may be caused, for 
example, by electromagnetic waves impinging on a non- 
linear medium or by mechanical vibrations. The elec- 
tromagnetic production of photon pairs by time-varying 
media and its connection with squeezing were analyzed 
by us in the general case some time ago [l|. A closely 
related problem is that of oscillating cavity walls (mir- 
rors). This problem, named the dynamical Casimir ef- 
fect, was investigated for the first time by Moore 0. A 
thorough review paper on the dynamical Casimir effect 
by Dodonov Q contains a long list of references. Later 
this line of research was extended to include all kinds of 
nonadiabatic transitions such as, for example, a proposal 
by Scully et al. [1] to generate photons by accelerating a 
neutral atom in a resonant cavity. 

Gravitational particle production due to changes of 
the metric with time was predicted a long time ago by 
Schrodinger [5| and later pursued by many others (cf., 
for example, [a, 0, [& [9|, LLO|). This effect has been con- 
sidered mostly in the cosmological context — the time- 
dependent metric was usually taken in the Friedman- 
Robertson-Walker form. In contrast, we consider here 
only local perturbations of the metric, avoiding global 
cosmological issues. 

In the case of gravitating orbiting bodies with total 
mass M and angular frequency of the orbital motion Q, 
there are two independent dimensional parameters: say 
GMQ/c 3 and Hc/GM 2 . Therefore it is not possible to 



estimate the order of magnitude of the effect by dimen- 
sional analysis (it could be any function of these two pa- 
rameters) and a detailed calculation is necessary. 

The starting point of our investigation is the set of 
time-dependent Maxwell equations, 

d t B(r, t) = - V x E{r, t), d t D(r, t) = Vx H(r, t), (1) 

which are valid in classical as well as in quantum electro- 
dynamics. For space-time-dependent e and fj, they can 
be written in the form 



d t B(r,t) + Vx 
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d t D{r, t) - V x Bi - V,t ^ = -Vx 0(r, t) 
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The right-hand sides in Eqs. ^ are responsible for the 
creation of photons. This process is best described with 
the help of the S operator (given as a time-ordered ex- 
ponential) , 



S = Texp(~ /d 3 rdtff/(r,£) 



(4) 



Hi(r,t) 
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D 2 (r,t) 



I3(r,t) 
2mo 



B 2 (r,t), (5) 



where D(r,t) and B(r,t) are the operators of the 
free field. In what follows we shall use the Riemann- 
Silberstein vector [III, 03 , 



F(r,t) 



D(r,t) , B(r,t) 



drk* 



fUv 



'Electronic address: birula@cft.edu.pl 



2e " V2Mo" J " ""V ( 2 ^) 3 
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where e(fe) = et(fe) is the normalized complex polariza- 
tion vector corresponding to the left-handed circular po- 
larization (positive helicity). The vector e*(k) = e^(fe) 
describes the right-handed polarization. The operator 
ajj(k) annihilates the left-handed photon with the wave 
vector k, and a R (k) creates the right-handed photon. 
The interaction Hamiltonian ([5} expressed in terms of F 
reads 

ff / = |(F + Ft) 2 -^(F-Ft) 2 . (7) 

Anticipating the smallness of the whole effect, we shall 
consider only the lowest order of perturbation theory: 



s -1--^ \ d 3 rdtH I (r,t). 



(8) 



II. RADIATION PRODUCED BY ORBITING 
DIELECTRIC OR MAGNETIC SPHERES 

As an introduction to the main subject, we consider the 
electromagnetic radiation produced by a homogeneous 
dielectric or magnetic sphere in a uniform motion on a 
circular orbit. In the case of a dielectric sphere of radius 
a, the function a(r,t) can be written in the form 



a(r,t) = ng(r- R(t)), 



(9) 



where g(r) = 9(a — r) refers to the sphere at rest, R(t) — 
R(cos(flt) , sm(fli) , 0) describes the orbital motion, k — 
(l/e r — 1), and e r is the relative static dielectric constant. 
For a magnetic sphere, the only difference is that in all 
the formulas e r is to be replaced by \x r . 

Assuming that initially there are no photons, the rele- 
vant terms of Hj are only those containing the products 
of two creation operators, i.e., 



S\0) « |0) 



4(2tt) 3 



J2 d 3 k! d 3 k 2 



Ai,A2 

x \/uJiLU2a(ki + k 2 ,coi +^ 2 ) 
xe^(fe 1 )-e^(fe 2 )4 i (fe 1 )4 2 (fe 2 )|0), (10) 

where a(k,u)) is the four-dimensional Fourier transform 
of a(r,t), 



&(k,uj) — / d re 
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g{r) j dte lult e lkR{t) 



a m (k,uj) 
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a m {k, Lo)2Ttd(uj — mil), 



K i m e m *f(ka)J m {k ± R), 



(11) 



(12) 



fej_ and <p are the polar coordinates of k in the orbital 
plane, k = \k\, and 

f(ka) — —j , 3 [sin(fca) — kacos(ka)] « 1 T7T~' (■"■^) 



The decomposition into Bessel functions in order to iden- 
tify the contributions with different frequencies has been 
used before by Nienhuis [l3| in his study of a rotating 
lens. 

According to formula (fT0|) , an orbiting dielectric sphere 
radiates pairs of photons. Note, that the Planck constant 
does not appear in Eq. (I10|) . even though the process 
of photon pair creation is a purely quantum-mechanical 
effect. This effect is not predicted by classical theory. 
The Planck constant will reappear in the formula for the 
radiated power. 

It follows from Eq. (T1TJ) that the sum of photon ener- 
gies is a multiple of Ml. Therefore, the sum over m is 
effectively restricted to positive values of m. The prob- 
ability w\ 1 \ 2 {k\ 1 k2) for emitting a pair of photons per 
unit time is 

,. . x 1 |e Al (fc 1 )-e A2 (fc 2 )| 2 
w\ 1 \ 2 (k 1 ,k 2 ) = y— j 



(2tt) 5 
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x ^wi^^a 6 f 2 (Ka)J^(Kj_R)5(w - mfl), (15) 



m— 1 



where K — k\ + k 2l tu — uj\ +lu 2 , and K± = kj_i + k± 2 . 
We employed here the standard trick 2ttS(ui — co') 2 — ► 
T5(lo — uj'), and we also used the relation 



|e Al (fc 1 )-e A2 (fe 2 )| 2 -(lTcos0) 2 /4. 



(16) 



The minus (plus) sign should be taken for the same (op- 
posite) polarizations and 9 is the angle between the vec- 
tors fci and k 2 . The presence of 1 =F cos^ shows that 
photons with the same polarization predominantly travel 
in opposite directions whereas photons with opposite po- 
larizations travel together. However, the probability of 
photons traveling in opposite directions is suppressed due 
to the small values of the Bessel functions near zero. 
Moreover, the factor lo\lo 2 reaches its maximum when 
u>i = lo 2 , since u>\ + lu 2 is fixed. Therefore, the most 
likely phenomenon is the emission of a photon pair with 
equal wave vectors and opposite circular polarizations. 

The sum over all Bessel functions represents contribu- 
tions from pairs of photons whose total energy is mKl. 
However, the intensities of higher harmonics are strongly 
suppressed. In order to quantify this statement, let 
us note that the arguments of the Bessel functions in 
Eq. (|15[) cannot exceed mRQ,/ c = mvu, where 



vr = Ril/c 



(17) 



is the orbital velocity measured relative to the speed of 
light. Since for small vr the ratio of two consecutive 
terms is approximately equal to 



(fca) 3 



10 



(m+l) 2 J^ +1 ((m+l)^) 
rrfiJ^mvn) 



1 



1 

m 



2m+2 



, (18) 



the contributions from higher harmonics can all be ne- 
glected and we can keep only the fii term. 

The angular distribution of the emitted pairs is de- 
scribed by the Bessel function J 2 (KR cos x), where x 1S 
the angle between the direction of the photon pair and 
the orbital plane. Since KR is small, the angular distri- 
bution is governed by cos 2 x- Therefore, the photons are 
emitted predominantly in the orbital plane. No photons 
are emitted in the direction normal to the plane. 

The main contribution to the total rate for emitting a 
pair of photons is 



1447T 3 



■I d 3 k 1 d 3 k 2 (l + cos 2 



x U1W2 f(Ka) Ji(Kj_R)S(w - Q). 



(19) 



To obtain an order-of-magnitude estimate of w, we intro- 
duce dimensionless photon wave vectors l\ = cki/Q and 
h = ck 2 /fl: 



w 



727 



^(|) v% d%d%(i + co S *e) 



x hl2f(Lv R a/R)Jl(L ± v R )S(h + l 2 - 1), (20) 

where T is the orbital period. For a nonrelativistic orbital 
velocity f 2 (Lv R a/R) « 1 and Jf(L±v R ) w (L^v R ) 2 /A. 
Hence, the total rate can be written in the form 



A 



k ( a\° a 

W ~288^t{r) Vi ' 

x / d 3 kd 3 l 2 (l + cos 2 e)hl2L\8{lx+l2-l). (21) 

The integral in second line is a pure number of the order 
of 1. Therefore, the total power of emitted radiation is 
approximately equal to 



Pf 



288vr 2 



e 8 m 

R) Vr ~T' 



(22) 



For a nonrelativistic motion, this value is extremely 
small, as was to be expected. It is worth noting that 
the probability of photon emission by a magnetic orbit- 
ing sphere differs only in having k 2 = (l/e r — l) 2 replaced 
by (l/p r - l) 2 . 

To compare the photon production by the orbiting ob- 
jects, caused exclusively by their electromagnetic proper- 
ties, with a similar effect caused by a time- varying gravi- 
tational field, we shall calculate now the radiation by two 
dielectric spheres orbiting around their center of mass. 
The radiation of this system due to the time-dependent 
gravitational field will be discussed in the next section. 
We assume that the bodies with masses Mi = pM and 
M 2 = (1 — \i)M and radii a\ and a 2 move on circular or- 
bits with frequency £1. Their trajectories are given then 
by the formulas 

Ri(t) = (1 - p,)R(cos{ilt),sin(Q,t),0), (23a) 

R 2 (t) = -^i?(cos(m),sin(m),0), (23b) 



where R = (GM/9, 2 ) 1 / 3 . In this case, 

a m (fe,w) = i m e mv M[^f{ka l )J m ((l - n)k x R) 
Pi 



K 2 



+ (-l) ro (l - n)—f(ka2)J m (ULk x R)] , 

pi 



(24) 



where p\ and p 2 are the average densities. The domi- 
nant contribution still comes from &i, unless the mate- 
rial coefficients are fine-tuned so that K\j p\ — K 2 /p 2 . In 
all the remaining cases, we may replace the difference 
Ki/pi — K2/P2 by some effective value R/p: 



|ai(fe,w)| 



2, ,2 »2 2h ,, ,_> / x - 

4" \r 



, _ M 2 kiR 2 p 2 (l-p) 2 R 



(25) 



The counterpart of formula (f2Tj) has the form 

A 



MV(l-/i) 2 (R\ 2 v* 



512tt 4 T \pj R 6 
x d 3 l 1 d 3 l 2 (l + cos 2 8)l 1 l 2 L 2 L 5{l 1 +l 2 -l). (26) 

Therefore, the radiated power is approximately equal to 
M 2 M 2 (l-/i) 2 (R\ 2 v R W. 



Pk 



512tt 4 



R e T 



(27) 



In the special case, when Ki/p\ = K 2 /p 2 = n/p, the 
leading term in the contribution from m — 1 vanishes 
and the contribution from m — 2 becomes dominant, 



\a 2 (k,uj' 



2 _ Af 2 fc 4 i?V 2 (l-M) 2 (^ ' 
* 64 



(28) 



and the radiated power, as compared to the generic case 
([2"7| . decreases roughly by a factor of v R , 



III. RADIATION PRODUCED BY A 
TIME- VARYING METRIC 

A double-star system creates a time-dependent gravi- 
tational field. We shall consider an isolated double-star 
system that modifies the metric only in its neighborhood. 
The electromagnetic radiation produced by such a system 
will be analyzed from the point of view of an observer lo- 
cated far away, where the space is flat. 

Of course, the exact form of the metric field for a 
two-body system is not known. However, for a ordinary 
stars (even for neutron stars) the Schwarzschild radius 
r is much smaller than the geometric radius and we 
can use the weak-field approximation. In this approxi- 
mation the diagonal components ho and h\ of the exact 
Schwarzschild metric in isotropic coordinates (cf., for ex- 
ample, [li|) 

ds 2 = h a c 2 dt 2 - h^dx 2 + dy 2 + dz 2 ), (29) 



ho = 



1 - r /4r 
1 + r /4r 



hi= 1 



Ar 



(30) 



can be replaced by 



To 

r 



/ii«l 



r 



(31) 



In this approximation, the theory is linearized and we can 
simply add the contributions to the metric from both 
bodies. For two orbiting stars, the metric components 
outside the stars are 



ho(r,t) = l- 
fti(r,i) = l + 



n 



l'2 



r-Ri(t)\ 
ri 



■*a(*)l' 

r-2 



(32) 

(33) 



|r-Jli(t)| |r-Jfc(t)|' 
where n and r 2 are the Schwarzschild radii of the stars: 



2GM 



n = fi- 



r 2 = (l-M) 



2GM 



(34) 



We have chosen a particular form of the metric tensor, 
but the calculated photon production rates are invariant 
under a change of the gravitational gauge. In the lin- 
earized version of gravity this change is given by (cf., for 
example, (la j) 



5g^{x) = dp£ u (x) + d v £n(x). 



(35) 



The proof of the invariance of the S operator under this 
transformation is given in the Appendix. 

The Maxwell equations derived from the variational 
principle based on the Lagrangian 



L = 



£o 
4* 



-gg^g^Uuhp, 



(36) 



written for an arbitrary metric, have the same form (fTj). 
but the constitutive relations are modified. In our case, 
for the diagonal metric, these relations are (cf., for ex- 
ample, [la ]) 



E(r,t) 



1_ h (r,t) 
«o V hi{r,t) 



D(r,t), 



H(r,t) = - 
Therefore, 

a{r,t) = /3{r,t) ■■ 



! h (r 7 t) 
hi(r,t) 



B(r,t). 



'h (r,t) 



-1 



V hi(r,t) 

r\ r 2 



\r-Ri(t)\ \r-R 2 (t)\ 
The interaction Hamiltonian is now equal to 
Hi =aF-F\ 
and the counterpart of formula (flT)]) is 



(37) 
(38) 

(39) 
(40) 

(41) 



5|0) « |0) 



d 3 fci 



d 3 k 2 



(2tt) 3 

X y/U>lU)2 &(kl + k 2 ,U>i +UJ 2 ) 

xel(k 1 )-e* R (k2)a{{k 1 )ai(k 2 )\0}. 



(42) 



Thus, the emitted photons always have opposite circular 
polarizations. The Fourier transform of a(r,i) is 



a(k,oj) = 
+ (l-/x) 



r>a± — oo 



rf 3 ? 



Ak-r 



dte iut e -ink-R(t)y 



(43) 



r>a 2 



Its decomposition into harmonics [cf. Eq. (|1 1[) ] gives 
16tt 2 



a m (fc,w) 



-i m e lm f> 



r\ cos(fcai) J m ((l — /J,)k±R) 



+ {-l) m r 2 cos{ka 2 )J m (t-ik ± R) 



(44) 



The probability of a pair production per unit time is 



w(fei,fe 2 ) = (2 \ 5 2 |ej?.( fc i)- e £(fe2)| 2 

OO 

x y^ |q m (iir, tju)| 2 (5(ct — mfi). 



(45) 



As in the previous case, the photons are predominantly 
emitted in the orbital plane with the same wave vectors. 
In the present case of a nonrelativistic motion of two 
orbiting bodies, due to a cancellation of the leading 
terms, \ai(K,w)\ 2 can be neglected as compared to 
\a 2 {2K \2w)\ 2 . The ratio 



\ai{K,w)[ 
\a 2 (2K,2zu) 



, ,'1 -2/j, 2 a{-ai K 



R 2 K ± 



(46) 



is small, of the order of v%, except for K± <C K, but this 
region contributes very little to the total rate. 

Keeping only the m = 2 term, we obtain the following 
final approximate formula for the total rate (after the 
substitutions k\ = 2£lli/c and k 2 = 2Vtl 2 /c) 



(»4 



w 



^(l-« 



7T 2 T 



x d 3 hd 3 l 2 (l + cost 



■ l\l 2 L ^ 



5(h + l 2 -l), (47) 



where we have used the relation (GMQ,/c 3 ) — v 3 ^. The 
integral, similar to the one in Eq. (|2ip . is a pure number 
of the order of 1. Therefore, the total power radiated by 
gravitating stars caused by the time-varying metric is 



64 2/1 , 2 10 2hn 

Pm ~ ^ (1 - fi) Vr —■ 



(48) 



For a double star of two solar masses and orbiting period 
of 1 h, we obtain vr — 0.0026 and therefore 



P M « 5.4 x 10" 27 



2htt 
~T~' 



(49) 



This effect is exceedingly small — one has to wait more 
than 10 22 yr for a pair of photons to be emitted. For 
solar-type stars the gravitational effect is weaker by three 
orders of magnitude than the effect due to the time- 
dependence of the dielectric constant. However, the radi- 
ation due to a disturbance of the metric might be stronger 
than the radiation due to the electric (or magnetic) prop- 
erties of the same double-star system. The ratio of the 
radiated power due to the gravitational effects ([48]) to 
that due to the electromagnetic effects ([2~T[) is 



M 



6.5 x 10' 



(I) 



2 G 8/3 M 2/3 



(50) 



This ratio varies quadratically with the density. When 
one of the stars is replaced by a typical neutron star with 
the density of about 10 17 kg/m 3 , the gravitational effects 
dominate over the electromagnetic effects by 25 orders of 
magnitude. 

It is also worth stressing that the orbital velocity ap- 
pears in Eq. (|4"8"[) to the tenth power. Therefore the effect 
would increase dramatically in the relativistic case. A 
significant increase of the velocity will change the power 
of emitted radiation by many orders of magnitude. 

In order to complete this investigation we shall com- 
pare it now with the power radiated by the same system 
in the form of gravitons. The emission of gravitons, like 
the emission of photons by an external current, is a purely 
classical phenomenon and we can expect that the associ- 
ated power Pq will be much larger. Indeed, the formula 
for Pq in the case of two gravitating bodies moving on 
circular orbits reads in our notation (cf., for example, 
0,0) 



64tt 7 Mc 2 
Pa = —v R — 



(51) 



Thus, the power radiated by double stars in the form of 
gravitons is enormous as compared to the power radiated 
in the form of photons due to the time-varying metric. 

In principle, at the more fundamental level, the pro- 
duction of photon pairs by gravitating bodies could be, in 
our opinion, described as a two-step process. In the first 
step, as predicted by linearized quantum gravity, the sys- 
tem produces a lot of gravitons. In the second step, col- 
liding gravitons produce pairs of photons (a single gravi- 
ton cannot decay) due to a coupling of gravitons to pho- 
tons [described by the interaction Hamiltonian (T4"l"j) ] . We 
believe that this mechanism will give results equivalent 
to those obtained by our semiclassical treatment of grav- 
itational effects. However, the calculation of the photon 



production rates along these lines would be much more 
complicated. 



APPENDIX A 

Our proof of the gravitational gauge invariance is pat- 
terned after a proof of the electromagnetic gauge invari- 
ance in quantum electrodynamics. In QED a change of 
the S operator under a change SA^ of the external elec- 
tromagnetic field is 



SS 



iT d 4 x [f{x)S]5A tl {x). 



(Al) 



This fundamental formula may even be viewed as a defi- 
nition of the four-current through the linear response of 
the system to a change of the external field. When the 
variation SA^, is only a change of gauge, 



SA^x) = d^A, 



(A2) 



the integral in Eq. (|A1|) . after integration by parts, van- 
ishes due to current conservation, 9 M j M = 0. 

In the case of gravity, the formula describing the 
change of S under a change of the metric tensor, a coun- 
terpart of Eq. (|A1[) . reads 



SS = — T 
2 



d 4 x[T^(x)S}Sg^(x), 



(A3) 



where T^ u is the energy-momentum tensor describing the 
coupling of matter to gravity. Inserting into this equation 
formula ([55)1 for the change of the metric tensor due to a 
change of the coordinates, we obtain after integration by 
parts 



5S = iT / d A x d„ [T>* v (x)S] ^(x). 



(A4) 



This integral vanishes because in linearized gravity 
d^T^" = 0. In our proof we argue that because the 
energy-momentum tensor is conserved, the S operator is 
a scalar, invariant under the transformations of the coor- 
dinate system. This is the reverse path to that followed 
often in general relativity (cf., for example, [181]). where it 
is argued that the energy-momentum tensor is conserved 
because the action is a scalar. 
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